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Abstract. We enlarge our set of supersymmetric models and update accelerator constraints in 
our precise calculation of the relic density of the lightest neutralino, which includes relativistic 
Boltzmann averaging, subthreshold and resonant annihilations, and coannihilation processes 
among charginos and neutralinos. 

1. Introduction 

The lightest neutralino is one of the most promising candidates for the dark matter in the 
Universe. A linear combination of the superpartners of the neutral gauge and Higgs bosons, 
it is believed to be the lightest stable supersymmetric particle in the Minimal Supersymmetric 
extension of the Standard Model (MSSM). 

In the near future, high precision measurements of the dark matter density may become 
possible from high resolution maps of the cosmic microwave background, and this may lead to 
constraints on supersymmetry. It is therefore of great interest to calculate the relic density of 
the lightest neutralino as accurately as possible. 

As a major step towards a complete and precise calculation valid for all neutralino masses 
and compositions, we include all concomitant annihilations (coannihilations) between neutralinos 
and charginos, properly treating thermal averaging in presence of thresholds and resonances in 
the annihilation cross sections (for details see [ |l| ) . 

2. Formalism 

Consider coannihilation of N supersymmetric particles with masses mi and statistical weights 
gi (first studied in ref. [ Normally, all heavy particles have time to decay into the lightest 
one, which we assume stable. Its final abundance is then simply described by the sum of 
the densities n = X) n «- When the scattering rate of supersymmetric particles off the thermal 
background is much faster than their annihilation rate, n obeys the evolution equation [ |j 
dn/dt = —3Hn — (o~ c gv)(n 2 — n 2 q ) with effective annihilation cross section (o~ c fiv) = A/n 2 q . 
The numerator A is the total annihilation rate per unit volume at temperature T, and n eq is 

the total equilibrium density. We find [ § A = (T/16vr 4 ) ds^Js - Am\ W(s) K^/T) and 
* Presented by P. Gondolo. 
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re eq = (T/27T 2 ) J2i 9i m iK^rrii/T). Here Ki(x) is the modified Bessel function of the second kind 
of order i, and W(s) = J2ij 9i9j&ij X(s,m 2 ,m 2 )/2^ X(s,m 2 ( ,m 2 c ) with X(x,y,z) = x 2 + y 2 + 
z 2 — 2xy — 2xz — 2yz. W(s) is a Lorentz invariant annihilation rate per unit volume in which 
coannihilations appear as thresholds at -y/s equal to the sum of the masses of the coannihilating 
particles. The independence of W(s) on temperature is a remarkable calculational advantage in 
presence of coannihilations: in fact it can be tabulated in advance, before taking the thermal 
average and solving the density evolution equation. 

The previous equations generalize the result of Gondolo and Gelmini [ [I] to coannihilations. 

3. Results 

To explore a significant fraction of the MSSM parameter space [ we keep the number of 
theoretical relations among the parameters to a minimum. We assume GUT relations for gaugino 
masses, keep only the top and bottom trilinear soft supersymmetry-breaking parameters, and 
use a single mass parameter for the diagonal entries in the sfermion mass matrices at the weak 
scale. We perform many different scans in parameter space, some general and some specialized 
to interesting regions. We keep only models that satisfy the experimental constraints on the Z° 
width, on the b — > 57 branching ratio, and on superpartner and Higgs boson masses. Here we 
show results for ~85,000 supersymmetric models that satisfy accelerator constraints as of March 
1998 (we include the ALEPH bounds at </s = 184 GeV [ @). 
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Figure 1. Neutralino relic density versus neu- 
tralino mass. Horizontal lines show the cosmo- 
logically interesting region 0.025 < £l x h 2 < 1. 



Figure 2. Neutralino masses m x and composi- 
tions Zg/(1 — Z g ) for cosmologically interesting 
models. 
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We obtain analytic expressions for the many Feynman diagrams contributing to the two-body 
cross sections at tree level for neutralino-neutralino, neutralino-chargino and chargino-chargino 
annihilation. Then for each set of model parameters, we sum over particle polarizations and over 
initial and final states numerically and tabulate the annihilation rate W(s). Thermal averaging 
and integration of the density equation finally give the neutralino relic density Q x h 2 = m x no/p CI it 
in units of the critical density p cr it- 

We find that coannihilation processes are important not only for light higgsino-like neutralinos, 
as pointed out before in approximate calculations [ Q], but also for heavy higgsinos and for mixed 
and gaugino-like neutralinos. Indeed, coannihilations should be included whenever \fi\ < 2\Mi\, 
independently of the neutralino composition. When ~ coannihilations can increase or 

decrease the relic density in and out of the cosmologically interesting region. 

Fig. |l| shows the neutralino relic density £l x h 2 versus the neutralino mass m x . To avoid 
artificial bands and holes in the point distribution, we have divided the plotted region into 
cells and marked those in which we find at least one model allowed by present accelerator 
constraints. The horizontal lines limit the cosmologically interesting region where the neutralino 
can constitute most of the dark matter in galaxies without violating the constraint on the age 
of the Universe. We take it to be 0.025 < 0, x h 2 < 1. 

Fig. |2] shows the cosmologically interesting region in the neutralino mass-composition plane 
{Z g is the gaugino fraction). This region is limited to the left by accelerator constraints, to the 
right by Q x h 2 < 1, and below and above by incompleteness in our survey of parameter space, 
except for the hole at 85-450 GeV where Sl x h 2 < 0.025. 

The main effect of coannihilations is to shift the higgsino region to higher masses. In particular 
the cosmological upper bound on the neutralino mass becomes 10 TeV. Differently from previ- 
ous approximate results [ , we find a cosmologically interesting window of light higgsino-like 
neutralinos with masses around 75 GeV. 

We conclude that if coannihilations are properly included, the neutralino is a good dark matter 
candidate whether it is light or heavy, and whether it is higgsino-like, mixed or gaugino-like. 
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